A new electromagnetic mode in graphene 
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A new, weakly damped, transverse electromagnetic mode is predicted in graphene. The mode 
frequency ui lies in the window 1.667 < < 2, where p is the chemical potential, and can be 

tuned from radiowaves to the infrared by changing the density of charge carriers through a gate 
voltage. 

In the past two years, a great deal of attention has been attracted by the discovery of graphene, a truly two- 
dimensional (2D) electronic system P, |3] (for recent reviews see Refs. @, 0]). Graphene is a monolayer of carbon 
atoms, and the band structure of electrons in graphene consists of six Dirac cones at the corners of the hexagon-shaped 
Brillouin zone @, Q, with the massless, linear electron/hole dispersion with the effective "velocity of light" V ~ 10 8 
cm/s. The special spectrum of the charge carriers leads to a number of interesting trans por t properties, which have 
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been intensively studied in the literature, see e.g. Refs. @, S, ©, M, ED, El El El El 
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Electrodynamic properties of graphene have been studied in Refs 
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Frequency-dependent conductivit y 
2ll |22| , thermoplasma polaritons [2 




11, 12, 18, a 



21, m 



i m las 129- 

and collective excitations of the graphene layer - plasmons 



have been theoretically investigated. It has been shown that the specific band 
structure of graphene leads to a certain modification of the plasmon and plasmon-polariton spectra, as compared to 
conventional 2D electron systems with the parabolic dispersion of electrons. 

In this Letter we show that the Dirac spectrum of electrons leads to a radically new feature of the electrodynamic 
response of the electron-hole plasma in graphene, as compared to conventional electron systems. We predict the 
existence of a transverse (TE) electromagnetic mode in graphene, the mode which cannot exist in systems with the 
parabolic electron dispersion. The new mode propagates along the graphene layer with the velocity close to the 
velocity of light, has a weak damping, and its frequency is tunable across a broad frequency range from radiowaves 
to the infrared. These properties may have a strong potential for future electronic and opto-electronic applications of 
graphene. 

To explain the essense of our finding, we briefly discuss the nature of electromagnetic modes in conventional elec- 
tronic systems. Isotropic and uniform three-dimensional plasmas can support (in zero magnetic field) both longitudinal 
and transverse electromagnetic modes. The electric field vector in the longitudinal (transverse) wave is parallel (per- 
pendicular) to the wave vector. In the 2D electron gas, however, only the longitudinal (or transverse magnetic, TM) 
modes - 2D plasmons and plasmon-polaritons - may exist under standard experimental conditions [271 ] . The spectrum 
of electromagnetic modes, propagating along and localized near the 2D electron gas layer, has the form 
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1mo(iS) y/ q 2 — uj 2 /c 2 







for the TM waves fegfl, and 



c 2 \J q 2 — lu 2 /c 2 







(1) 



(2) 



for the transverse electric (TE) waves [29j]; here cr(u>) is the local dynamic conductivity of the 2D gas and c is the 
velocity of light. As seen from ([T]) and @, the TM (TE) modes may exist if and only if the imaginary part of <r(u>) 
is positive (negative). In conventional 2D electron systems, realized, for instance, in GaAs/AlGaAs quantum-well 
structures, the conductivity can be described by the Drude model cr(u>) — in s e 2 /m(uj + ry), where n s , e, m and 7 are 
the density, the electric charge, the effective mass and the scattering rate of 2D electrons, respectively. As ct"{lo) > 0, 
only the TM waves (plasmons, plasmon-polaritons) can propagate in such structures. 

In graphene with the massless Dirac form of the electron/hole dispersion the situations is different. To demonstrate 
this, we begin with the Hamiltonian H = Va a p a , which determines the energy spectrum of charge carriers in the 
electron (I — 2) and hole (I = 1) bands, 

Eu = (-l) l hVk, 1 = 1,2, (3) 

and the corresponding wave functions |kZ). Here a = (x.y), a a are Pauli matrixes, p a is the momentum operator, 
and k = (k x ,k y ), k = ^Jk 2 + k 2 is the vave vector. Using the self-consistent-field approach 30] (alternatively, the 
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FIG. 1: (Color online) The dynamic conductivity of the graphene layer, in units e 2 g s g v /16h, as a function of the frequency 
= huj/fj, in the collisionless limit at zero temperature T/fi — 0. Imaginary and real parts of the intra-band and inter-band 
contributions, as well as their sum (the total conductivity) are shown. 



Kubo formalism [8(, or the random-phase approximation [31|) we calculate the high-frequency conductivity of the 
system <j a p(uj) = a{ui)5 a 0. Scattering is ignored (i.e. 7 = 0) here. The conductivity consists of the intra-band and 
the inter-band contributions [12|, LL8| , 
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where /(-Ekz) is the Fermi distribution function, v a — Va a is the velocity operator, and S is the sample area. In the 
gate controlled graphene systems, where the density of charge carriers can be tuned by the gate voltage, we obtain at 
T/ fx — > 0, that the intra-band conductivity Q assumes the Drude-like form, 

in s e 2 V 2 _ e 2 g s g v 4i 
amtra{Uj) ~ (a, + i0)/i " "W^' (6) 



while the inter-band contribution gives 



Winter (w) = — *- - (|fi| - 2) In 
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see Figure [TJ Here Q = ?uj/li, /1 is the chemical potential, n s — g s g v ^i 2 / AttT% 2 V 2 is the density of electrons at T = 0, 
and g s = 2 and g„ = 2 are the spin and valley degeneracies in graphene (six Dirac cones give the valley degeneracy 
g v — 2 because only one third of each cone belongs to the first Brillouin zone @). 

As seen from Eqs. (|6]) and (|7|), as well as from Figure [IJ the imaginary part of cri ntr a{^) is positive, like in the 
standard 2D electron systems, so that the intra-band contribution alone could not provide the conditions of the TE 
mode existence. The imaginary part of the inter-band contribution, however, is negative and diverges logarithmically 
at — > 2. This divergency is associated with the step-like behavior of the real part a' inter (ui), which describes the 
inter-band absorption of radiation at tku > 2/i. Using ((6|) and ([?])■ the TE-mode dispersion |2]) can now be written as 



^^ = M 1 ±(n 2 + a \ 



(8) 
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FIG. 2: The value of Q 2 — Q 2 as a function of the dimensionless frequency fi for the TE mode. 



where Q = hcq/ ' \x. Notice that apart from the dimensionless wave- vector Q and the frequency f2, Eq. ([8]) depends 
only on the fine structure constant a = e 2 /Tic. The TE mode © doesn't decay at T — and exists in the window 
1.667 < f2 < 2, where the term in brackets in ([5]) is positive. To give a numerical example, at n s ~ 6 x 10 12 cm -2 the 
TE mode should exist in the window 115 THz < / < 140 THz, at n s ~ 10 11 cm~ 2 - in the window 15 THz < / < 18 
THz, and so on. Choosing an appropriate value of the gate voltage, one can tune the TE mode to a desired frequency 
range. This may be useful in designing devices for infrared, terahertz, and microwave electronics. The TM mode (2D 
plasmon-polariton) does not exist at 1.667 < ft < 2; its dispersion relation (Q} has real solutions only at SI < 1.667. 

Figure illustrates the dependence of \J Q 2 — SI 2 on the frequency $1, Eq. ([5]). Due to the small factor a — e 2 /tic in 
the right-hand side of ([8]), the deviation of the wavevector Q from the frequency SI in the wave is small as compared 
to the values of Q and £1 themselves. This means that the TE mode propagates along the 2D graphene layer (x 
direction) with the velocity close to the velocity of light, u) < cq, and that the localization length of the wave in the 
perpendicular (z) direction is much larger than its wavelength in the x direction. 

The frequency dependence of the conductivity at finite temperature is shown in Figure [3J At T > 0, the real part 
cr'(w) becomes finite at Q < 2, and the TE mode acquires a finite damping, Q = Q' + iQ" . This damping, however, 
is very small. Using that a<Clwe get from Eq. ([8]) 



where a is the normalized conductivity of the layer, a = (e 2 g s g v /16fr)cr. Figure [4] exhibits the factor a'(us)[— a"(ui)] 
as a function of O at different temperatures. One sees that a'(uj)[— a"(cu)} is smaller than 1 and hence that q" /q' does 
not exceed 10 -4 even at T ~ 0.1/i. At the density n s ~ 6 x 10 12 cm~ 2 the Fermi energy [i corresponds to T ~ 3000 
K in graphene, so the TE mode should be easily observable at room temperature. 

To conclude, we have predicted a new transverse electric mode in graphene. The existence of the mode is directly 
related to the linear "relativistic" spectrum of charge carriers in graphene; in a conventional system of 2D electrons 
with the parabolic dispersion such a mode cannot exist. The mode frequency is widely tunable across the range from 
radiowaves to infrared frequencies, depending on the density of electrons or holes in the system. The damping of the 
mode is very weak even at room temperatures. 

The work was partly supported by the Swedish Research Council and INTAS. 
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FIG. 3: (Color online) The real and imaginary parts of the dynamic conductivity of graphene, in units e 2 g s g v /16ft, as a function 
of frequency Q — Tioj/ fj, at finite temperatures. 
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